Strongly nonlinear magnetization above Tc in Bi2Sr2CaCu208+5. 
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Using high-resolution magnetometry we have investigated in detail the magnetization M above 
the critical temperature Tc in Bi2Sr2CaCu208+5. In a broad range of temperature T above Tc, we 
find that M{T, H) is strongly non-linear in the field H. We show that as T — s- Tc, the susceptibility 
x{T,H) diverges to very large values (x ~* "1) if measured in weak H. In addition, M{H) displays 
an anomalous non-analytic form AI ~ H^^^ in weak fields with a strongly T-dependent exponent 
S{T). These features strongly support the proposal that, above Tc, the pair condensate survives to 
support significant London rigidity. 

PACS numbers: 74.72.-h,74.25.Ha,74.25.Dw 
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I. INTRODUCTION 



The notion that superconductivity in the cuprates is 
destroyed by thermally created vortices has gained sub- 
stantial experimental support. In single crystals, evi- 
dence for vortex excitations that persist above the criti- 
cal temperature Tc to a temperature Tnnset ~ 120 K has 
been obtained from the Nernst effect An enhanced 
diamagnetic signal that scales accurately as the Nernst 
signal has also been observed using high-field magnetom- 
etry 0. These results reveal that, above Tc, significant 
pair-condensate amplitude survives up to intense mag- 
netic fields H. Further, in ultra-thin films, the kinetic 
inductance has been observed to persist above Tc jsj . The 
sensitivity of superconductivity to phase fluctuations in 
cuprates has been investigated theoretically by several 
groups 0,0. 



In low- Tc superconductors, the condensate vanishes 
above Tc, except inside evanescent droplets created by 
amplitude (or Gaussian) fluctuations. The diamagnetic 
susceptibility x from Gaussian fluctuations is very small 
(x ^ 10~^) y. Given the vortex scenario in cuprates, the 
magnetization M measured above Tc ought to be qualita- 
tively different from the Gaussian picture. However, re- 
solving the signal in tiny crystals is very challenging, and 
the available results are confusing 0JH0)0|- An exper- 
iment performed on aligned grains Q claims agreement 
with the mean-field Gaussian theory. Other experiments 
are seemingly consistent with the vortex scenario, but un- 
cover anomalies that are hard to interpret . Here we re- 
port detailed measurements in Bi2Sr2CaCu208 (Bi 2212) 
showing that the weak-field susceptibility x diverges ex- 
ponentially consistent with Kosterlitz-Thouless (KT) be- 
havior. The strong nonlinearity observed in x suggests 
the existence of a subtle rigidity in the pair wave function 
^ which persists above Tc. 
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FIG. 1: Magnetization AI vs. H measured in Sample 1 in 
fields \H\ < 1,000 Oe (Panel a) and in fields \H\ < 40 Oe (b). 
M is measured in the SQUID magnetometer in the long-time 
averaging mode (~12 h per curve). In Panel (a), the sharp 
notch near H = QiorT <Tc {= 86 K) is Hd. These features 
are shown more clearly in Panel (b) with a 25-fold expansion 
in the field scale. The dotted line is AI/H = x/(l + ^cx) = 
2.55 with X = -0.95 and Nc = 0.64 for Sample I with H||c. 
The peaks locate Hd. 



II. EXPERIMENTAL DETAILS 



To cover a broad range of fields, we have combined 
torque magnetometry with SQUID magnetometry. In the 
torque experiment, the crystal is glued at the tip of the 
Si cantilever beam with its c-axis at a tilt-angle if ~ 15° 
to H. In our experiment, the torque r is the sum of 
a paramagnetic term from the spin moment nip and a 
diamagnetic term from the magnetization M of interest 
here, viz. f — [nip -I- VM] x iiqH, with V the sample 
volume and /io the permitivity. For temperature T > Tc, 
M is strictly along the z-axis (we take z||c). An in-plane 
com pon ent Mx is resolvable only below the irreversibility 
line [ll| (from hereon, Mz = M). We define the effective 
magnetization Mcff = t / {^iqHxV) where Hx = Hsinip. 
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FIG. 2: Curves of M vs H in Bi 2212 (Sample I) measured by 
torque magnetometry for T above Tc = 86 K. The pronounced 
curvature evident in Fig. ^ persists to '^lOO K in fields upto 
2 T and higher. The maximum beam deflection angle Stp is 
0.15°, or lO'V- 

For (/5 <g; 1, we have 

M,ff{T, H) - Axpi/. + MiT, H,), (1) 

where the anisotropy = Xc — Xa is the difference 

between the c-axis and in-plane spin susceptibihties. As 
explained in Ref. , the very weak T dependence of Axp 
in cuprates allows the strongly T-dependent M{T, H) to 
be extracted reliably from Me//. All curves reported 
here are fully reversible (hystereses caused by pinning 
of vortices are only seen below ~50 K). We studied 2 
slightly underdoped Bi 2212 crystals with Tc = 86 and 
85 K (Samples I and II) with very similar results. 

If M is nearly linear in H (valid if T > 1.1 Tc and 
H > 0.5 T), Eq. n implies that the torque varies si- 
nusoidally, viz. t/V — /io(Axp — x)HxHz = Aisin2(p, 
where Ai{T,H) ~ and changes sign above Tc, as ob- 
served in cuprates (see below). 



III. MAGNETIZATION AND SUSCEPTIBILITY 

First, we discuss the curves of M{H) below Tc- The 
curves of M{H) in Fig. ^ were measured by SQUID 
magnetometry. Much analysis have focussed on the curve 
of M{H) at the "crossing temperature" Tg (= 84 K) at 
which M{H) is very nearly H independent. By detecting 
the flux-exclusion in very weak fields, we have determined 
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FIG. 3: The field dependence of M{T,H) in Bi 2212 {Tc = 
86 K) from = 10 Oe to 20 T at T = 79-115 K (Sample I). 
The log-log plot shows that M{H) obeys Eq. Has J/ -» 0. 
The slight upturn at 83 and 84 K near 10 Oe reflects Hci (Fig. 
□a). The plot combines SQUID results (10 to 1000 Oe) and 
toque magnetometry results (500 Oe to 2 T). Torque results 
up to 20 T are also shown at selected T. The dashed line (of 
slope 1) is the linear response of Eq. I^at 87 K. 



that the 3D transition at Tc (86 K) actually lies 2 K 
above Tg. In the interval T^-Tc, a prominent anomaly of 
the curves of M{H) is that they curve steeply towards 
zero as 77 — > even though full flux exclusion obtains 
when H < Hd (the lower critical field). This feature 
may be seen in earlier reports on M 12], but has not 
received comment. Closer inspection reveals that M{H) 
displays a notch feature which we show in expanded scale 
in Fig. ^p. The notch corresponds to the initial entry 
of vortices at Hd- Below Hd, flux exclusion gives a 
linear M-H variation with a steep slope given by M/H = 
x/(l + Ncx), where Nc ~ 0.64 is the demagnetization 
factor. Taking x — —0.95, we have M/H = —2.55 which 
is plotted as the broken line in Fig. Both the full 
flux-exclusion implied by % ~-l and the sharp onset are 
evidence for a high degree of electronic homogeneity and 
uniformity of the condensate. 

In low-Tc superconductors, \M\ falls steeply above Hd, 
and continues its monotonic decrease until H attains the 
upper critical field Hc2- Here, the behavior is qualita- 
tively different. Just above Hd, \M{H)\ falls to a min- 
imum, but then slowly increases over several decades in 
H. Following the initial entry of vortices, the sample 
steadily becomes more diamagnetic with increasing B. 
As we show, this paradoxical pattern reflects the phase- 
disordering nature of the transition at Tc in hole-doped 
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FIG. 4: The temperature and field dependence of the sus- 
ceptibility x{T,H) = M{T,H)/H. In weak H, \x\ increases 
by over 5 decades from 100 K to Tc, reaching the full flux- 
expulsion value (enhanced by demagnezation) shown by the 
arrow. The dashed curve is a fit to Eq. |31with b = 1, d — 
12 A, a = 10 A, and Tkt ~ 84 K. Open circles represent the 
KT correlation length ^kt obtained from the fit. 



cuprates. 

We turn next to temperatures above where the 
resolution of the torque cantilever 10~^ emu at 10 T) 
allows M to be studied in detail. Figure El displays the H 
dependence of M at temperatures increasing from near 
Tc to 100 K. Over this broad interval, M increases with H 
with a curvature that becomes increasingly pronounced 
near Tc- These curves are strikingly similar to those in 
the interval Ts-Tc, except for the conspicuous absence of 
the notch feature associated with a finite Hd ■ 

We may combine the SQUID and torque results over 
our full field and temperature range in Fig. |3| The 
SQUID results cover the \ow-H regime (10 Oe to 1000 
Oe) while the torque curves extend to much higher fields. 
In this combined plot, M initially increases as a power 
law in H, reaches a broad maximum at fields 1-10 T and 
then decreases monotonically. The fan-like dispersion of 
the straight lines in low fields implies that M has the 
power-law behavior 

M{T, H) = A{T)H^/^^'^^ [H ^0), (2) 

described by a strongly T-dependent exponent 6{T) (A 
is field- independent). In the interval 84-105 K, where 
S > I, the non-analytic form of M implies that, as 
H ^ 0, the diamagnetic susceptibility increasingly ex- 
ceeds linear-response behavior (dashed line drawn for 87 
K). 



FIG. 5: The T dependence of the weak-field magnetization 
exponent 5{T) in Samples I (solid circles) and II (open circles). 
The reciprocal S{T)~^ is plotted as solid and open triangles 
for I and II, respectively. As T ^ , 5{T)~^ decreases 
smoothly to 0. In the shaded region from Ta to 105 K, where 
5 > 1, linear magnetic response is absent even at 10 Oe (Eq. 



In Fig. H the susceptibility x{T,H) = M(T,H)/H 
is displayed versus T, with H as the parameter. Evi- 
dently, x{T, H) is highly sensitive to H especially near 
Tc- If X is measured at a fixed H, say 0.5 T, it stays 
non-divergent across Tc. However, this is not intrinsic. 
Instead, from 105 to 86 K, |x| measured in a very weak 
field increases by 5 decades to attain values 0.1-1 just 
above Tc- As discussed later, this steep increase reflects 
a rapidly changing phase correlation length. 

Previous experiments on "fluctuation diamagnetism" 
in cuprates failed to observe a divergent x because the H 
applied was too large The strong H dependence 

in x{Ti H) now call their conclusions into question. Fur- 
ther, vortex pinning effects, especially strong in polycrys- 
talline samples of La2_xSrxCu04 and YBa2Cu307_5, in- 
troduce extrinsic features in M which render analyses 
problematical [l^. 



IV. KT CORRELATION LENGTH 

There are 2 major aspects of the non-linearity in 
M (T, H) above Tc- The first reflects the divergent growth 
of the phase correlation length ^(T) while the second is 
the singular behavior in Eq. |2] The 2D phase correlation 
length measures the region within which phase coherence 
prevails in each layer. In KT theory long-range phase co- 
herence is destroyed by the appearance of (anti) vortices 
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at the KT transition (at Tkt)- Hence the KT corre- 
lation length S^iiT is the average spacing between these 
thermally created vortices. Above TKT^ M{T, H) is lin- 
ear in H in low fields. The weak-field susceptibility is 0| 

XKT = -^^eKT. {T>Tkt.H^QI (3) 

where fc^ is Boltzmann's constant, d the bilayer spacing, 
and 00 the superconducting flux quantum. The KT cor- 
relation length is exponential in the reduced temperature 
t' = T/Tkt ~ 1, viz. = ae^/^, with 6 1 a non- 
universal constant and a a cut-off length ~ the vortex 
core size. 

Clearly, an interlayer coupling needs to be included 
to describe the actual 3D transition at Tc (see below). 
Nonetheless, the KT theory describes the essential fea- 
tures of a divergent £,{T) and the roles of thermally- 
induced vs. field-induced vortices in each layer. In Fig. 
01 the dashed curve is a fit of Eq. El with 6=1.0 and a — 
10 A. The fit captures quantitatively the broad features 
of the 5-decade rise in xi^i H) in weak TJ, matching its 
curvature both near Tc and near 100 K. 

The KT correlation length ^kt inferred from the fit in- 
creases from '^8 nm at 105 K to 700 nm at 86 K. A second 
length scale is the spacing gb — yj4>o/ B of vortices in- 
serted by H. To observe the intrinsic divergence in S^kTi 
we need as > £,kt or B < Bth, where Bth = (^o/^kt 
is a field-scale set by the density of thermally created 
(anti)vortices. As T ^ in a fixed field B <C Bth, x{T) 
initially increases rapidly (and Bth decreases). However, 
X saturates to a constant when Bth approaches B, as 
evident for the curves in Fig. 01 

V. SINGULAR RESPONSE AS if ^ 

Despite the reasonable fit to x{T), Eq. O does not 
account for the non-analytic behavior expressed in Eq. 
121 As mentioned, the deviation between M{H) at 87 K 
and the dashed line in Fig. grows as H decreases. We 
now discuss in more detail the T dependence of this non- 
analyticity. From 105 K to 84 K, the weak-field slopes 
in the log-log plots in Fig. decrease smoothly from ^^1 
to zero, i.e. S{T) diverges to values exceeding 20. The 
steepness of this increase is evident in the plots in Fig. 
m for Samples I and II (solid and open circles, respec- 
tively). The continuous nature of the increase in 5{T) 
is made clearer by plotting its reciprocal, which falls be- 
low 1 and smoothly approaches as T ^ T+ (triangles) . 
In Fig. O the shading highlights the interval from Tg 
to ^^105 K in which 6 exceeds 1. In the shaded interval 
above Tc, x diverges without limit in the limit H — >0 
assuming Eq. [3 remains valid. Hence, despite the ab- 
sence of a true Meissner state, a feeble H induces a very 
large diamagnetic screening current. Below T^ however, 
this limit is interrupted at Hd by the appearance of the 
Meissner state, as discussed in Fig. 



Generally, the Meissner state originates from the rigid- 
ity of the pair wave function In a field, this "London 
rigidity" suppresses the paramagnetic current response, 
leaving only the full diamagnetic current [l9|. Here, our 
magnetization results imply that, despite the loss of the 
Meissner effect, ^ retains some rigidity above Tc. How- 
ever, the rigidity is very fragile and observable only in 
a feeble field. The existence of a singular diamagnetic 
response above Tc has not been predicted (for e.g., S is 
strictly 1 above Tkt in KT theory). Singular behavior 
in the weak-field vortex-liquid response is discussed in 
Ref. 113, and in Ref. 



VI. DISCUSSION 

We now return to the anomalous behavior of M{H) 
in the interval Tg-Tc where full flux-exclusion prevails for 
H < Hci- As mentioned, the curves below Tc closely 
resemble those above (compare Figs. ^ andjH). Aside 
from fields below Hd, the curves of M in the log- log 
plot in Fig. O assume a fan-like dispersion that smoothly 
extends the behavior in this interval to above Tc- The ex- 
ponent 1/J(T), which measures the slopes, continuously 
decreases across Tc (Fig. EJ. This continuity strongly 
implies that the anomalous magnetization in the 2-K in- 
terval below Tc and above Tc share the same origin; both 
represent the diamagnetic response of a strongly phase- 
disordered 2D superconducting state with a very high 
depairing field iJc2 > 80 T 0. Between and Tc, the 
magnetization is comprised of a weak-field 3D Meissner 
state that pre-empts the phase-ordering transition of the 
2D condensate at T^. Application oi H > Hd destroys 
the Meissner state and uncovers the underlying 2D state 
whose magnetization continues to grow with field. This 
accounts for the unusual "minimum" in \M\ just above 
Hd mentioned in Fig. The steep fall of \M\ just 

above Hd reflects the destruction of the 3D Meissner 
state, while the subsequent power-law increase reflects 
the behavior of M in the underlying 2D condensate which 
is robust to intense flelds. 

A hint of our results may be seen in previous torque 
experiments |l[ll|. In crystals of TlzBazCuOg+a with 
Tc = 25 and 15 K, the measured curve of r vs. Lp displays 
a puzzling, unexpected term Ai sin 2ip which is actually 
dominant above Tc [TH . Using Eq. ^and our analysis, we 
now understand this term to be the robust 2D magnetiza- 
tion which persists to intense H and high T as discussed 
above. As noted, the rapid decrease in \M\ above Tc, 
compared with the mild change in A^p, engenders a sign 
change in Ai{T,H) as seen (Fig. 5 in Ref. [ll|). 

The transition in cuprates has been frequently com- 
pared to the 3DXY transition with large anisotropy a^^, 
where a = J^/ J with J the intralayer coupling and Jj^ 
the interlayer coupling 0, 0, 0. Our findings that 
the 2DKT theory gives a broadly accurate description of 
X(T, H) and that Tc-Tkt ~ 2K < confirm that Bi 
2212 falls in the extreme limit a <C 1. [By way of com- 
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parison, in the planar ferromagnet K2CUF4 with a ^ 1, 
the KT theory also describes well M vs. H above its 
Curie temperature Tc 0|. However, the 3D transition 
Tc = 6.25 K pre-empts the 2D KT transition {Tkt ~ 5.6 
K) JJJ. Below Tc, the destruction of the 3D Meissner 
state here obviously has no parallel in the planar magnet.] 
Unfortunately, the limit a ^ 1 is intractable analytically. 
Currently, numerical simulations of M on lattices ,1_5, 16] 
are not accurate enough to compare with our measure- 
ments, but we hope the detailed measurements here will 
motivate simulations on larger lattices. 

As we cool from 300 K, the first signs of the pair con- 
densate gradually appear at Tonset ~ 120 K as a vortex- 
Nernst signal and an enhanced diamagnetic suscepti- 
bility ^5] . The phase of ^ is strongly disordered by vortex 
motion. However, the continued growth of these signals 
to intense fields (30-45 T) implies that the depairing field 



Hc2 hes much higher (80-100 T). Below -105 K, the non- 
linear behavior of x becomes increasingly apparent. In 
weak iJ, x{T, H) undergoes a 5-decade increase as Tc is 
approached, whose broad features can be understood in 
terms of a diverging phase-correlation length ^, as pre- 
scribed in KT theory (Fig. 0)). In the shaded interval 
in Fig. El X is further enhanced by a singular power-law 
dependence suggestive of increasing London rigidity. 
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